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A closed expression is obtained for the cross-section for Coulomb excitation of 
levels of the giant dipole resonance of given angular momentum and phonon number. 
Applications are made to the Goldhaber- Teller and Steinwedel- Jensen descriptions 
of the resonance, at non-relativistic and relativistic bombarding energies 
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The giant dipole resonance (GDR) is one of the best-studied collective modes of nuclear 



excitation. 



to protons 



t has been modeled, macroscopically, as a bulk oscillation of neutrons relative 
l|, or as local isovector fluctuations of neutron and proton fluids It can 
also be modeled, microscopically, in terms of isovector linear combinations of particle-hole 
excitations of the nuclear ground state (jjj, 4, 5, 6]). It has been studied experimentally using 
reactions induced by gamma rays, light ions, and by the sharp pulses of electromagnetic 
radiation associated with projectile nuclei moving at relativistic speeds, i.e. relativistic 
Coulomb excitation. A broad survey of the GDR and other resonances, and different different 
excitation methods, is given in ref.p]]. In this paper, we will consider relativistic Coulomb 
excitation of GDR states which are described by macroscopic models. 

In many situations, the de Broglie wave length for the relative motion of the projectile 
and target in a Coulomb excitation experiment is small compared to the linear dimensions 
that characterize the system. Then a semi-classical approach can be used, in which this 
relative motion is described in terms of a classical orbit, whereas the internal changes of 
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the projectile and target are described using quantum mechanics. We are interested in a 
situation in which the projectile remains unexcited, but the target is excited to the states 
associated with the GDR. 

A useful first approximation to an oscillatory situation, such as the GDR, is to assume 
that the restoring forces are proportional to the displacement from equilibrium. With this 
approximation, the GDR is dynamically equivalent to an isotropic 3-dimensional harmonic 
oscillator. This assumption leads to the familiar harmonic oscillator spectrum, in which 
eigenstates are characterized by phonon number, total angular momentum, and angular- 
momentum-z-component. If the oscillator picture were exact, all the (N + 1)(N + 2)/2 
eigenstates with N phonons would be degenerate in energy. Deviations from the oscillator 
picture would lift this degeneracy, but if the deviations were spherically symmetric we would 
still have the 21 + 1-fold degeneracy of the angular-momentum eigenstates with I = N, N — 
2, AT- 4,..., or 1. 

The coupling between the electromagnetic pulse due to the projectile and the internal 
oscillating degrees of freedom associated with the GDR of the target can usefully be ap- 
proximated by an expression that is linear in these oscillating degrees of freedom. If this 
approximation is made, an exact solution can be found for the Schrodinger equation that 
describes the time evolution of the target jfj, Formulae have been published in the liter- 
ature for the total excitation probability of all GDR states of given phonon number, when 
the relative motion of the target and projectile is along a specified orbit. In this paper, we 
decompose this total excitation probability into the contributions of phonon states of given 
total angular momentum. For example, we show how to find the excitation probabilities of 
four-phonon states of angular momentum 0, or 2, or 4, whereas the previously published 
formula yielded only the total four-phonon excitation probability. 

The model studied in this paper is highly simplified, since it uses a pure oscillator de- 
scription of the GDR, and a linear approximation of the coupling between the projectile 
motion and the GDR degrees of freedom. More realistic calculations have been done, for 
e Xam ple in ref, gflflfiul which have proven nsefn! in the ana lysis of a wide variety of 
Coulomb excitation experiments. We believe that our simplified model may, nevertheless, 
be useful in indicating some general trends, especially with respect to the variation of ex- 
citation cross-section with bombarding energy, and the way this depends upon the angular 
momentum of the final state. 
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Unfortunately, these GDR phonon states of specified angular momentum are not clearly 
resolved in the excitation spectra. Indeed, superposed on the multiphonon GDR states are 
collective excitations of other characters, such as giant quadrupole and giant octupole exci- 
tations (see, e.g., ref. {3]). Thus we cannot check our predictions for excitation cross-sections 
of GDR states of given angular momentum against any currently available data. However, it 
is possible that future measurements of angular distributions of the decay products of GDR 
states will give information about the angular momenta of these states. For example, the 
gamma rays emitted by the i = member of the two-phonon sextuplet will have a spher- 
ically symmetric angular distribution, whereas the five £ = 2 members will emit gamma 
rays with quadrupole and hexadecapole distributions. In situations such as this, it will be 
important to be able to predict the excitation cross-sections of N— phonon states of specified 
angular momentum. 



II. EXCITATION PROBABILITIES 

A. The coupled Equations 

The time-dependent Schrodinger equation for the perturbed target wave function is 

^^- = [Ho + V(t)m) (2.1) 

with H the unperturbed target Hamiltonian, whose eigenfunctions are ip a and eigenvalues 
e a . The solution of this equation is expressed in terms of occupation amplitudes a a (t) which 
occur in the expansion 

*(t)=J]a a (t)e-^Va. (2.2) 

a 

Substitution of (2.2) into (2.1) yields a set of coupled differential equations 

iha a (t) = Y^e^-'f* < ^ a \V{t)\^ p > a p {t). (2.3) 

P 

The intial conditions appropriate to a typical nuclear reaction are a a (— 00) = 5 a0 , corre- 
sponding to the requirement that the target be in its ground state at the start of the process. 
The probability that the reaction leaves the target in the final state ip a is then |a a (oo)| 2 . 
In a Coulomb excitation reaction, the perturbation matrix elements required in (2.3) are 
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HQ 

<^\V(t)\^ Q > = J [tpf(r,t) P p a (r) - ±Af(r,t)-J Pa (r)]d 3 r . (2.4) 

Here </^ et (r, i) and A^ et (r,t) are, respectively, the scalar and vector potentials associated 
with the electromagnetic field created by the charged projectile. The properties of the 
target states ip a and ipp are expressed in (2.4) by the transition charge density pp a {v) and 
current density J i g a (r). 

B. Excitation probabilities in a Cartesian basis 

In this paper we will be concerned with situations in which the target Hamiltonian H is 
that of an isotropic three-dimensional harmonic oscillator with reduced mass M and natural 
frequency u. The state labels used in Sec. II. A can be taken to represent the triplet of 
quantum numbers (n x ,n y ,n z ), specifying the numbers of oscillator quanta in the x, y and z 
directions. Furthermore, we will be working in the regime in which the interaction matrix 
elements (2.4) can be approximated by 

<Vv,n,,jV(t)|VH,n' s X> = -|^C, % ,jR)[F(i)'R+G(i).P]i,„,„,(R) (2.5) 

where R(= X, Y, Z) represents the degrees of freedom undergoing harmonic oscillations and 
P represents the conjugate momenta. 

Since both Hq and V separate in Cartesian coordinates, the problem reduces to one- 
dimensional forced oscillations in the X, Y, Z directions, and the occupation amplitudes 

a-n X) n y) nX t ) factorize 

a nx ,n y ,nM = KM K(t) &„,(*)■ (2-6) 
In particular, the b nx (t) satisfy 

ihb nx (t) = _^ e *(«--<)* < n x \F x (t)X + G x (t)P x \n' x > b K (t) 

+ e~** (]f^F x (t) - i][^G x (t)\ v^TT b n+1 (t), (2.7) 
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which must be solved with the intial condition b nx {— oo) = S nxt Q. It is not difficult to verify 
that (2.7) is satisfied by 



with 



a x (t) 



dt' 



w +j ,/^. (0 



&(*) = i / dt' 

-oo 



F x (t') nr^ 



— i 



2h 

G x (t') 



a x (t')e 



-wtf 



(2.8a) 

(2.8b) 
(2.8c) 



V2MUu V 2h 

from which it follows that (3 x (t) + f3*(t) = —\a x (t)\ 2 . Here and in the following, we shall 
use the notation on, without explicit time- dependence (t), to denote the quantities ai (oo); 
the same convention will be also used for the quantities ai and hi. Then the probability of 
populating the final target state ip nx ,n y ,n z is 



V 

I n x ,n v ,n z 



\Kf\h 



, 2 I6 I 2 

n y I I u n z | 
2\n z 



(i^i 2 r-(i^i 2 ^(i« 2 i 2 )"^ _ ( | 0j2+ i a „ |2+laz|2) 



Tl x .Tly .Tl z . 



(2.9) 



in which a y and a z are defined as in (2.8b), but using F y (t'), G y (t') and F z (t'), G z (t'), instead 
of F x (t'),G x (t'). Note that the'Toisson distribution" result (2.9) involves "on-shell" Fourier 
transforms of F x (t),G x (t), as given by (2.8b). 



C. Excitation probabilities in a spherical basis 

The Cartesian result (2.9) is well known [8j. However specifying the eigenstates of Hq in 
terms of n x ,n y ,n z is not as convenient as using principal and angular momentum quantum 
numbers n, £, m. The advantage of using n, £, m is that a spherically-symmetric deviation 
from a perfect harmonic oscillator Hamiltonian will not mix the states labeled by different 
(£, m), nor will it split the 2£+ 1 states with different m- values and the same £. However, the 
quantum numbers only useful for a perfect oscillator. Since we cannot expect 

perfection in a harmonic description of the GDR, but we can expect spherical symmetry, 
it would be advantageous to use oscillator eigenstates characterized by n, £, m rather than 

Tl'X) H"yi ^z- 

Using (2.2), (2.6) and (2.8) we can write the exact time-dependent target wave function 
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as 



i?^?J i?^2 



(2.10) 

where \n x ,n y ,n z > represents a harmonic oscillator eigenstate with n x ,n y ,n z quanta in the 
x, y, z directions, respectively. It is convenient to represent the oscillator eigenstates in terms 
of boson creation operators c+, c+, c+ acting on a normalized vacuum (ground state), |0 >: 

\n x ,n y ,n z > = / ^ — 10 > (2.11) 

All the properties of these oscillator eigenstates can be obtained algebraically, starting from 

< 0|0 >= 1 (2.12a) 

K, 4} = <W, [c M , c„] = [c+ c+] = (2.12b) 
c M |0 >=< 0|c+ = (2.12c) 
If we substitute \n x ,n y ,n z > from eq.(2.11) into ty(t) given by eq.(2.10), we get 

m = ^W+AW+AW-iH £ [^^^W4]"1e-^a,(t)4]^[e-^ 2 (t)4]^ |q ^ 

5 ?^y 5?^2 

p /3 x (t)+l3y(t)+l3z(t)-i%u>t ST^ [f ^ ( a z CO C z + + Q *(^) c z )] N | n 

AT! I U> 

jv=o 

= e l3 x {t)+Py{t)+p z {t)-i%^t x e e- 1 - t K(t)c++a H (t) C ++ az (t)c+) (2.13) 

We now use the polynomial identity 

6 = (2n)!!(2n + 2£+l)!! ^ ^«»^ b ) ( 2 ' 14 ) 

which is proven in Appendix A. Here the solid harmonic 3^m( a ) = a ^m^P) * s a homogeneous 
polynomial of degree t in a x ,a y , a z with the same rotational transformation properties as 
the spherical harmonic Y^(a). With the help of eq.(2.14), we can rewrite eq.(2.13) in the 
form 



\&(t) = e feW+A/ffl+fe(«)-'§^ e -iu>(2n+e)t 



x 



-l) n v / 4^ L 4==^==^=L y m 4 aJt), a z (t)) \n£m > (2.15a) 
' v/(2n)!!(2n + 2£+l)!! ml W ' yW * W> " Jl ; 
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where 

\( c +)2 i ( c +)2 , ( c +)2]n 

\ nim> = (_i)V4^L +UJ + l *J L y e m (c+,c+ 4)[0 > (2.15b) 

Since |n£m > is a homogeneous polynomial of degree 2n + £ in c+, c+, c+ acting on the 
oscillator ground state, it must be an eigenstate of the harmonic oscillator Hamiltonian 
with 2n + £ quanta. Its rotational transformation properties are the same as the spherical 
harmonic Y^. Moreover, it is shown in Appendix B that the numerical factor in (2.15b) 
guarantees that \n£m > is normalized. Therefore, eq.(2.15a) implies that 

< n£m\^(t) > = e /3*{t)+Mt)+f3z(t)-4ut x (_-gm 

r-[<*l(t) + a 2 At) + a 2 Jt)] n , 
x -1 VI^4 ^ vKJ z ULy e _ m (a x (t),a y (t),a z (t)) (2.16) 
V ; v/(2n)!!(2n + 2£+l)!! mV W wW ' W; V ; 

is the amplitude that \n£m > is occupied at time t. 

The quantity of interest for the interpretation of experimental data is the sum of the 
excitation probabilities of all states of given n, £: 

e 



2 + |cU 2 P 



(2n)!!(2n + 2£+ 1)!! 
x e-d-l 2+ ^l 2+ l-l 2 ) £ \y m (a x ,a y ,a z )\ 2 (2.17) 

In evaluating this sum, care must be taken when complex conjugating the solid harmonics 
because the arguments a x ,a y , a z may be complex. We have 

i 

m=-£ 

e 

= (-l) m £ ^m(«x. «D> «*)^4K> (2-18) 
m=-e 

For any two vectors ri and r 2 , the spherical harmonic addition theorem gives 

£ £ 
m=—£ m=—£ 

= [(sj + + + ^ + g±i x P, ( -= mi + ± »» ) (2.19) 
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The argument of the Legendre polynomial in (2.19) is just f 1 • f 2 . Equation (2.19) can be 
regarded as a polynomial identity in the six variables Xi,yi, z±, x 2 , y 2 , z 2 , which we can apply 
to the particular choice 

xi = a* x , yi = a*, zi = a* z 
x 2 = «!, y 2 = a y , z 2 = a z 

Then (2.18) and (2.19) yield 

1 

\yL( a x,®y,®z)\ 2 

m=~£ 

= |K) 2 + K) 2 + (a 2 ) 2 rx^P,( 
and eq.(2.17) gives 

T) 2£ + 1 , 2 2 , 2 \2n+£ 

Vn ' = (2n)!!(2n + 2f + l)!! l K) +M +M ' 

e-d-l 2 +l^l 2 +l-l 2 ) x Pl ( , ' a j 2+ J a C +KI ! 2l l (2-20) 







2 + 




2 + 




to 


IK) 2 + ( 


«,) 2 + 


(a 2 


) 2 I 



for the total excitation probability of states with given n, £ 

As a simple example of the application of (2.20), consider the two-phonon states {n = 
0,£ = 2) and ( n =l,£ = 0). Using P (x) = 1, P 2 (x) = (3x 2 - l)/2, we find that 

Pi,o = i|K) 2 + (%) 2 + (a,) 2 | 2 e-d-l 2+ l^l 2+ l-l 2 ) (2.21a) 

(|a x | 2 + K| 2 + K| 2 ) 2 - ||K) 2 + K) 2 + (a 2 ) 2 | 2 
' 0,2 - r 



x e -(KI 2 +KI 2 +M 2 ) (2.21b) 

The sum of (2.21a) and (2.21b), the total excitation probability for two-phonon states, is 
seen to be 

<p 10 + <p 02 = (Kl 2 + \ a v\ 2 + \ a z\) x e -(M 2 +KI 2 +M 2 ) 

2 

which is the same as the sum of V nx>nytnz of (2.9) for the six n x , n y , n z combinations for which 

n x + n y + n z = 2 

In the application of (2.20) to Coulomb excitation of the GDR, we will find that a x = 0, 
a y is real and a z is pure imaginary. This is because the influence of the projectile in the 
z-direction is opposite at ±t, whereas its influence in the y-direction is the same at ±t. This 
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different time behavior leads to different behavior under complex conjugation of the Fourier 
transforms that determine a y ,a z . The significance of the opposite signs of (a y ) 2 and (a z ) 2 
is evident in (2.20). 

In general, a x , a y and a z will be functions of the impact parameter, b, which characterizes 
the projectile orbit. Therefore V n ,i will also be a function of b. The excitation cross-section 
involves an integral over impact parameter, 

POO 

<r n>i = 2ix V n/ (b) b db. (2.22) 

" t>min 

The lower limit, 6 m i n , is of the order of the sum of the radii of the projectile and target 
nuclei [yj]. Because the electromagnetic pulse due to the projectile becomes more adiabatic 
as b increases, V n ^{b) decreases strongly for large b, and the upper limit of the integral in 
Equation (2.22) can be safely taken to be of the order of a few hundred Fermi. 



III. APPLICATIONS 



A. Application to the Goldhaber- Teller model of the giant dipole excitation 

The Goldhaber- Teller |l| model is based upon a division of the target degrees of freedom 
into three sets: 

p[, p' z locating the protons relative to the proton mass-center, 

n.j_, n'jy locating the neutrons relative to the neutron mass-center, 

R locating the proton mass-center relative to the neutron mass-center. 

States of the giant resonance excitation are postulated to have the form 

Vtf = X (p , i,...,Pz;n / 1 ,...,n' w )$^(R). (3.1) 

All states have a common spherically symmetric intrinsic state x(p'd ■■■■P'z'i n 'n •••• n 7v)' but 
they differ in the relative motion of the neutron and proton mass centers, as specified by 
the different functions $^ £ (R), which describe small oscillations with an approximately 
harmonic restoring potential. 

The transition charge and current densities between two states and are 

Pnim,n'i'm'( r ) = e J d Pl, dp Z ] d^, dn N 

z 

x ^ 5(r - Pi )*™ £ (pi, p z ; ni, n N )V*f (p x , p z , n u np$.2a) 

i=i 
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Jnem,n'e' m '(^) = 7. : / dp±, dpz', diii, dn N 



2m p i 



z 

x 

i=i 



Zj 

<H r - Pi) ( ^™ £ (Pi> Pz; n 1; n JV )V Pi *^'(pi, p z ; n 1? njv ) 
7 (pi, ...pz; ni, n jV )V Pi ^(pi, Pz ; n 1; njv )) , (3.2b) 



where m p represents the proton mass. The vectors pi, ...pz', n i, n N locate the tar- 
get protons and neutrons relative to the target mass center. They are related to the 
p[, ....p' z ; n^, ....n' N defined previously by 

p, = fR + p' t (3.3a) 

Uj = (/-l)R + n;., (3.3b) 

where / = -=^-= (see Figure 1). 

If we use (3.3) and (3.1) in (3.2a), we get 



Pnem,n'e'm'(r) = e J dRdp' v ....dp' z , dni, ....dri N \xk(p'i, -, Pz', n 'i> •••> n ' 

z 

*;lR)^'(R)E^ r -/ R -ft) 



z 

x $™ £ (R)^f'(R)' 

i=i 

The structure of this expression suggests we define an intrinsic density function 

z 

Pint(s) = e / rfp / 1 ,....rfp / z ;rfn , 1 ,....rfn / 7V |x(p / 1 ,....p / z ;n / 1 ,....n / iV )| 2 ^5(s-p^), (3.4) 
^ i=i 

in terms of which p n £ m ,n , e'm'(*') and Jn£ m , n 'f m '(r) can be written 

Pnem,n'e>m>(r) = J d 3 R Pint (r - /R)$^(R)<r (R) (3.5a) 

W«W = / ^Pint(r-/R)(^(R)V R ^'(R)-^'(R)V R ^(R)) 

(3.5b) 

The physical significance of Pmt(s) is seen from Figure 1 to be the charge density defined 
relative to the mass center of the proton distribution. 

The amplitude of the oscillation in R is determined by the parameter v = Mcu/h, where 
M is the reduced mass associated with the relative oscillation of the proton and neutron 
mass centers, and Tiuo is the characteristic energy of the oscillation. For the giant dipole 
resonance, the oscillation energy is approximately 

kuj ~ 79 A' 1/3 MeV 
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(Q). Thus 

NZ uj NZ „ huj NZ 79A- 1 / 3 MeV NZ „ _ 2 

z/ = — — m p - = —m p ^— ~ — — x 939 MeV x — — ~ 1.91— 7- fm 

A h A (he) 2 A (197.3 MeV fm) 2 A 4 /3 

The amplitude of the oscillation of the proton mass center is then 

/X 7^ = f X 7^"^yf fm - (3 - 6) 

If this distance is small compared to the distance over which Pmt(r) changes by an appreciable 
fraction of itself, such as the thickness of the nuclear surface, then it should be useful to do 
a Taylor expansion of Pmt(r — /R) about R = 0: 

p int (r-/R)~p int (r)-/R-V r p int (r) + ••• (3.7) 

This enables us to approximate (3.5a and b) by 

Pnem,n>t>m>(r) - J d 3 R[p int (r) - fR • V rPint (r)]$;f (R)<f (R) 

= 5n,n'V WPmt(r) - /V rPint (r) ■ J $™ £ (R)R<f(R)d 3 # (3.8a) 

Wvm'tt * I d 3 it!p in t(r)[^(R)V R <f(R)-<f(R)VR^(R)] 

Pint(r) 



Zm p 



J ^ J f(R)P^ n J(R)d 3 R (3.8b) 
If we combine these expressions with (2.4), we get 

<^\V(t)\^ n J > = 5 n , n ,5 u ,5 m , m , J d 3 r^\r,t)p int (r) 

- / J rfV c ct (r,t)V rPint (r)- J ^(R)R^'(R)d 3 R 

- ^ J d 3 rAf(r,t)p int (r) • J $^(R)P<f (R)d 3 R (3.9) 

The first term is a monopole integral, independent of the internal degrees of freedom of the 
target. Its effect can be absorbed into a time- dependent phase factor multiplying the wave 
function, and we ignore it in the following discussion. We then see that the remaining terms 
in (3.9) are of the form (2.5), with 

N f 

F(t) = - J dV^f(r,i)V r p int (r) (3.10a) 

G(t) = — ^- / d 3 rA^(r,t)p int (v)z . (3.10b) 
ZrripC J 
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To estimate the validity of the approximation (3.7), let us consider the particular example 
of a 40 Ca target. Then Equation (3.6) yields f/y/u ~ 0.724/40 1/3 fm ~ 0.212 fm for the 
amplitude of the oscillation of the proton mass center. Since the proton charge density 
is approximately constant from its center out to the surface region, whose thickness is 
approximately 1 fm, we see that the amplitude of the GDR oscillation is indeed small 
compared to the distance over which the charge density changes by an appreciable fraction 
of itself. Thus we would expect that Equation (3.7) is a reasonable first step in the analysis 
of GDR data. 

We have not used the term "electric dipole approximation" to label (3.7), since that term 
usually implies a comparison of a wavelength with the size of a charge distribution. This is 
not the comparison we need to justify (3.7). 

We now distinguish between regimes in which the projectile is relativistic or non- 
relativistic: 



1. Relativistic projectiles 



In this approximation Uj, the projectile linear momentum is very large compared to the 
transverse impulse the projectile receives as it moves past the target. Then the trajectory of 
the projectile can be approximated by a straight line, along which the projectile moves with 
constant speed v. We choose our axes with the target center at the origin, the projectile 
trajectory in the y — z plane, a constant distance b from the z axis, and the projectile moving 
in the z direction. Then ^^(r, t), A ret (r, t) are the Lienard-Wiechert potentials 121 ] 

•yZpe 



x 2 + (y — b) 2 + 7 2 (^ — vt) 2 



A ret (r,t) = -^ et (r,t)z 

c ° 

We can now use (3.10) and (2.8b) to calculate the a x ,a y ,a z needed in (2.20): 

F(f) 



(3.11a) 
(3.11b) 



OL 



alt' 



+ ^ 



Mu 

~2h 



G(0 



Suit' 



h N 

2MwI 



d 3 V c et (r,^)V r p int (r) + 



where 



alt' 



Muoh v 
2 ZrripC 2 



dV c ct (r,a;)p int (r)f3.12) 
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is the "on-shell" Fourier transform of the retarded potential (3.11a). A convenient multipole 

nn 

expansion of this function is [111, [13J 



<pTm = ^E e "' mi ^ (^) £ iu (t'') i »( f '- < 313a > 

m \ I / x=\m\ ^ ' 



where the coefficients Qx, m are defined by 

■\-\-m / \ A— m 

^-^(pj (^)VM2A + l)(A-m)!(A + m)!X; 



(27) 2n (n + m)!n!(A - m - 2n)\ 
(3.13b) 

We assume that the target proton charge distribution is spherically symmetric, so Pi n t( r ) = 
p int (r). Then 

J nn9d6d<pY£(e, <p)Vp int (r) = 5 XA ^-p[ nt (r) x q(m) (3.14a) 
where q(l) = -^±M g(0) = z and q(-l) = and 



sin 6d6dipY* (6, (p)p in t(r) = 5 x ,o5 m ,o V An p int (r) (3.14b) 

If (3.11) is substituted into (3.10), and (3.14) is used to evaluate the angular integrals, the 
result is 

a x = (3.15a) 
IQNuj Z p e 2 fub\ f°° . /u , 2 



a = -itt^/ — — — / jo(-r) p int (r)r dr (3.15b) 

w m p ZAh 7?r \jv J Jo \c 



a z = it* 



32Nuj Z p e 2 T ^ fub\ 



rripZAh 7 2 t> 2 
Use has been made of the relation 



K ° (t^) / J ° (~c T ^ Pint ^ r2dr (3.15c) 



' h ( - r J Pint( r ) r2 dr = --/ jo ( -r) pmt(r)r 2 dr 

This result is gauge invariant, since it involves only on-shell Fourier components of the 
interaction. 



2. Non-relativistic projectiles 

Here Newtonian mechanics is used to obtain a hyperbolic trajectory for the projectile 
jloT ] . Following the conventions used in the previous Section, we choose the coordinate axes 
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so that the target center is at the origin, the trajectory is in the y — z plane, reflection- 
symmetric across the y axis, with the projectile moving in the direction of increasing z. If 
r'(t) locates the projectile center at time t, then the scalar and vector potentials are 

Zpe 



(p(r,t) 



|r - r'(t)| 



A(r,t) = ^(r,t) 



(3.16a) 
(3.16b) 



Here v(t) = dr'/dt is the projectile velocity at time t. 

The hyperbolic trajectory is characterized by the lengths of its semi-transverse and semi- 
conjugate axes, a and b, related to the asymptotic kinetic energy E and the angular mo- 
mentum £ by 



ZpZ^e 
2E 



2mE 



(impact parameter). 



(3.17) 



The eccentricity of the trajectory is e = + (b/a) 2 . The adiabaticity parameter is £ = 
aw/v, with v the asymptotic relative speed (v = y/2E/m). Then F(t)and G(t) of eqn.(3.10) 
are replaced by 



F(t) 



£ f^r Zpe 

A T J r 



r'(t)\ 



Pintl r ) r 



N T 2 (x x'(t)+yy'(t)+zz'(t)) 
-M ZfZt€ 



G(t) 



Zpe 



(3.18a) 



Z^r/ipC 2 
Zpe 2 



d 3 rp- mt (r) 



x x'(t) + y y'(t) + z z'(t) 



r ; (t) 



m p c* 



x x\t) + y ^(t) + z i'(t) 
r'(t) 



(3.18b) 



The time integrals corresponding; to eqn. (3.12) can be performed using the methods given 
by Alder and Winther [lOj. Some details are given in Appendix D. The result is 



a. 



CM ,, 



iZ Z — / 2N ^ huJ e -?f 

he y ZtAt^pC 2 



7 



7 2 - 1 



h-h 



Zp Zt - — a 



2N T hu 



he y Z-pATmpC 2 



ti b 
e ? 2 _ 



7 



2 r 

-<3 



7 2 - 1 e 



(3.19a) 
(3.19b) 

(3.19c) 
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where 



h 

h 
h 
h 



-f e cosh(ui) 



-£e cosh(ui) 



-f e sinh(ui) 



cosh(u>) 



cos(^w) + esinh(w) sin(£u>) 
1 + (e sinh(w)) 2 



dw 



cos((,w)dw 



sinh(w) 



— sin(£u>) + e sinh(w) cos(£w;) 



dw. 



(3.20a) 
(3.20b) 
(3.20c) 
(3.20d) 



1 + (e sinh(w7)) 2 

These integrals are easily calculated numerically. Note that in this non-relativistic approxi- 
mation, only the total target charge is relevant, not the radial charge density distribution. 



3. Overlap between relativistic and non-relativistic regions. 

It is of some interest to consider the transition between the regions of applicability of 
the relativistic and non-relativistic approaches to Coulomb excitation. The principal ap- 
proximation limiting the validity of the relativistic approach, as the bombarding energy 
decreases, is the assumption that the projectile moves along a straight-line trajectory at 
constant speed. Actually, at low speeds the projectile is deflected by the Coulomb field of 
the target, with a classical scattering angle of 9 = 2ArcCot (|) = 2ArcCot ( Zp 2 ^ g2 j ■ As an 
example, consider 208 Pb projectiles Coulomb deflected by a 40 Ca target. Then the deflection 
angle would be less than 1° if Eh > 82x20 2 xL44 x Cot(.5°) = 1.35 x 10 5 MeV fm. For a 
Coulomb-dominated trajectory, b > 12 fm, so we will not encounter trajectory deflections 
more than 1° if E > 1.13 x 10 4 MeV, or a projectile energy of 1.13 x 10 4 /208 ~ 54 MeV 
per nucleon. Moreover, for an assumed straight-line trajectory with h = 12 fm, the decrease 
in projectile speed as the Pb projectile moves past the Ca target is less than 1%. Thus we 
would not expect the deviations from a constant velocity orbit to cause significant problems 
for the relativistic approximation in the 208 Pb- 40 Ca system as long as the projectile kinetic 
energy per nucleon exceeds 54 MeV. 

Another aspect of the relativistic approximation is the neglect of the recoil of the target as 



the projectile moves past it. Winther and Alder ll|] have shown that this may be corrected 



to some extent by replacing the relativistic adiabaticity parameter ^ in eqn (3.15) by 
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This correction is applied in the results to be shown below. 

The most serious error introduced when the non-relativistic approximation is applied to 
fast projectiles is the omission of the relativistic sharpening of the electromagnetic pulse 
experienced by the target. It is this sharpening that causes the adiabaticity parameter to 
be rather than The presence of the I/7 factor increases the amplitudes of high- 
frequency components in the pulse, which are necessary for the population of a high-energy 
excitation mode such as the GDR. Thus we might expect the non-relativistic Coulomb 
excitation formalism to underestimate the GDR excitation cross-section as the projectile 
kinetic energy increases into the relativistic region. 

Figure 2 shows a comparison of the results of relativistic and non-relativistic calculations 
for the excitation of the one-phonon GDR level in 40 Ca, as a result of 208 Pb-induced Coulomb 
excitation. The curve labelled "relativistic" is certainly applicable at the high-bombarding- 
energy end of the energy range. The above qualitative considerations indicate that it may 
also be valid down to bombarding energies per nucleon of about 50 MeV. This inference is 
supported by the fact that at 50 MeV per nucleon, it agrees with the non-relativistic result, 
which should be valid at this lower energy. Thus there is good reason to believe that, in this 
case, the relativistic result is valid down to about 50 MeV per nucleon. At this energy the 
magnitude of the cross-section is so low that detection of the GDR would be difficult. Thus 
the relativistic theory seems to be valid over the entire useful energy range. 

Figure 2 also shows that the difference between the predictions of the relativistic and 
non-relativistic theories, as the bombarding energy increases, is in the expected direction. 



B. Application to the Steinwedel-Jensen model of the giant dipole excitation 



This is a model |2j in which protons and neutrons oscillate relative to each other, not in 
the bulk relative motion of the Goldhaber- Teller model but in local isovector fluctuations. 
Our version of this model follows the presentation of Greiner and Maruhn 15|. We present 
only the version of the theory in which the projectile is relativistic, because we have seen 
in the last section that this theory is applicable over the entire energy range for which 
multiphonon GDR levels are observable. 

Let s N (r, t) and s p (r, t) be neutron and proton displacement fields, i.e. when a fluctuation 
occurs a neutron that was in equilibrium at r moves to r + s JV (r, t). Because of the isovector 
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nature of the fluctuation, s N (r, t) and s p (r, t) can be expressed in terms of the single vector 
field s(r, t): 



Sp(r,t) = -js(r,t) 

The corresponding velocity fields are the time derivatives: 

. , ds p (r,t) Nds(r,t) N . . 

v r i = — = ±-L-L = — v fr t) 

A ' ' dt A dt ~ A 



v„(r,t) = 



ds„(r,t) 
9* 



Zds(r,t) 
A dt 



z , \ 
^v(r,t) 



(3.21a) 
(3.21b) 



(3.22a) 
(3.22b) 



Let n be the equilibrium nucleon number density, so that -jn and -^n are the equilib- 
rium proton and neutron number densities. As a result of the displacements (3.21), these 
number densities become 



Z Z 

n P (r,t) = -n (l-V-s p (r,t)) = -n + rj(r,t) 

N N 

n N (r,t) = — ra (l- V-s„(r,*)) = —n - r](r,t) 



(3.23a) 
(3.23b) 



where the isovector density fluctuation rj(r : t) is related to the isovector displacement field 
s(r, t) by 



r](r,t) 



NZ 



n V • s(r,i) 



(3.24) 



For small fluctuations, the kinetic and potential energy densities are given by 



^(r,0 = -K(r,t)(v p (r,t)) 2 + n N (r,t)(v N (r,t)) 2 } 



— n 



m p ZN 



ds x (r,t) 
dt 



ds y (r,i) 
dt 



ds z (r,t) 
dt 



(3.25a) 



U(r,t) = 4^(r,t)) 2 = 4^ , 

n n \ A z 



SZ\ 2 fds x (v,t) | ds v (r,t) | 9s M (r,t) y 



dx dy dz 

Here a s (~ 23 MeV) is the symmetry energy parameter. If the Lagrange equations of motion 



d dC 



d dC 



d dC 



d dC 



dt \ <9ff J dx\ d^f J dy \ d^f J dz \ d^f 



dx 



dz 



dC 

ds n 



are applied to the Lagrangian density C = T — U, the result is 



d 2 s 



dv 



n o m p-Q^ = n o™ p — = - 8a s Vi] 



(3.26) 



18 



Number conservation for the protons and neutrons 
drip . . Z drift 



-V • (npvp 



Z ^ 
-— n V • v P , 



dt A dt 

can be expressed with the help of (3.22) and (3.23) as 

dr] NZ 

m = -^ V ' V 

Combining (3.26) and (3.27) leads to the d'Alambertian equation 

d 2 r] 



-V • (n N v N ) = -— n V ■ v N 



(3.27) 



dt 2 



with 



u 



u 2 V 2 ?7 



I8a s NZ 
m p A 2 



(3.28a) 



(3.28b) 



The Steinwedel- Jensen model of giant dipole excitations [2j is based on solutions of (3.28a) 
of the form 

NZ hi^r) 



r](r,t) 



-n - 



-r • a(i) 



(3.29a) 



A 2 r 

where the time-dependent amplitude a(t) exhibits harmonic oscillations with frequency uj: 

d 2 &(t) 



dt 2 



- LU 2 &(t) 



(3.29b) 



To show that (3.29) satisfies (3.28), we observe that ^-a(t) is a linear combination of spherical 
harmonics Y^(r), and that 

2 



V 2 ji (^r) Y,l(f) 



(;)*(r)* 



The irrotational displacement field consistent with (3.29) and (3.24) is 

,2 



s(r,t) 



u- 



— - — r ■ a(t) 



(3.30) 



The boundary condition that determines u is the requirement that the radial component 
of s(r, t) should vanish at the nuclear surface. If this were not true, fluctuations would imply 
a complete separation of protons from neutrons on either side of the nuclear surface which 
is regarded as unphysical (although it occurs in the Goldhaber- Teller model jl]). Thus it is 
required that 



u 



= f ■ s(r,t)\ r=R = — 



d (uj \ ( uj \ d u . (uj \ 

~dr iu r ) r ■ a « + » L R ) Tr {T ■ aW) = ^ in ") T ' a(i) 



d 
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It is therefore required that be a zero of j[ (^-R). The lowest zero is ~ 2.08, and 00 is to 
be determined by the eigenvalue condition 

uj = 2.08^ (3.31) 
R 

So far the picture is one of classical harmonic oscillations of a two-fluid system within 
a closed spherical volume of radius R. To quantize this picture, we quantize the "virtual 
oscillator" (3.29b), i.e. we treat the coordinates of a quantum harmonic oscillator 

of frequency u. The states of the giant dipole excitation would be eigenstates of this virtual 
oscillator. To get the size parameter associated with this virtual oscillator we identify its 
potential energy ^Muj 2 (a x + a 2 +al) with the volume integral of the potential energy density 
(3.25b) 



/ d 3 rU(r,t) = f d 3 r{ V (r,t)) 2 = 4a s n 

Jr<R n Jr<R 



J r <R \ 

1 



A 2 



d 3 r ( J1Ku r ■ a 

r<R 



= -Muj 2 (a 2 x + a 2 + a 2 z ) 



A simple calculation shows that the size parameter of the virtual oscillator is 



% / 3hu* hu(2.08Ay 



Mu \l32na s n (^) 2 u3f 2 - \ qjl ( q )) 2 d q " V 8a, x 0.462' 



(3.32) 



Now the charge and current densities can be expressed in terms of the operators a x ,a y , a z 
and the conjugate momenta P x , P y , P z : 

p(r) = en P (r) = — n + erj(r) 

Z NZ jigr) 

= e ~r n o + e ~ nr n ° — ~ — r ■ a (3.33a) 

, , , , ZN , , ZN d , N 
J(r) = en P (r)v P (r) ~ e— n v(r) = e— n ^s(r, t) 

ZN u 2 f 3l ^r) \ ZN u 2 -fj^r) \ 

If these charge and current densities are substituted into (2.4), we get an expression of the 
form (2.5) (with R of (2.5) replaced by a here), from which we can extract 



F(t) = ~^n J dWc(r,t) J -^r (3.34a) 
G(t) = ^^o(^) 2 / tfrtf* (r,t)V (^r) (3.34b) 
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and 



a 



oo 

/ / dt' 

-oo 



F(t>) 



V2Mh 



CO 



Suit 



—I 



h eNZ 
2Muo A 2 



n / d 3 r<p T c et (r,uj) 



Mhu v eZN it 2 /\ 3 rct . 9 /Oi(^) 



With the help of (3.13), we find 



OL, 



.8rrNZ Z p e 2 n 



h K ( ^ 

2Mu K jv' 



x 



A 2 h^v 

r (irji(-r)ji(-r) 



u 



u f 2 , ■ \ ■ r u \ 
- I r drj 2 (—r) ] 2 {-r) 

v J u c 



(3.35a) 



(3.35b) 



a 7 



8nNZ Z p e 2 n 

"^42 hv 

if. 



2M^ K ^ X 



1 uv 
3 1 2 

c I" 2 1 ■ / w \ ■ /U \ 2u / 1 
-/ r drji(-r)ji(-r) +-- 1 + — 
v J u c 3 v \ 2j z 



,0J s ,cj 



jo(-r)j (-r)r dr 



u 



,w , . ,00 



R 

r 2 dr j 2 {-r)j 2 (—r) 
« c 



(3.35c) 



IV. HIGH BOMBARDING-ENERGY LIMIT 



According to Equation (3.15), the high-bombarding-energy behavior of a y and a z in the 
Goldhaber- Teller model [1] is determined by 



7 ^oo . , 32A Z p e 2 

OLn, > — ITT \ — *-r- 

Xl m p ZATiu cb 



7^00 
a z — ► 7T 4 



32 Ac; Z p e 2 / 7 c 
Ln — 



m p ZAh 7 2 6 2 \ub 



jo (~r) p- mt {r)r 2 dr (4.1) 
jo (^V) Pint(r)r 2 dr (4.2) 



Thus a y has a finite high-bombarding-energy limit, whereas a z approaches zero as Ln(7)/7 2 . 
The corresponding expressions for the Steinwedel- Jensen model |2j are obtained from Equa- 
tions (3.35): 

2 



a, 



.SttNZ Z p e 2 c 2 u 2 n / % 



A 2 hu 5 bR 



2Muj 31 \ r 



uR 



uR ( ujR 

4 cos I + 

u \ u 



ujR 
u 



4 sin 



uR 



11 



(4.3a) 
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8irNZ Z p e 2 c 2 no / h ( ub\ . / uR 



a * —* nT~ V 2M^ K ° [To) 31 1 — 



uR 2 , fuR\ 1 u 4 

Ji 



uR fuR\ ( fuR\ 2 \ fuR 



w -v«; 1 - - •■ - (4 - 3b) 

The eigenvalue condition that determines u is the vanishing of j[ (— )■ Therefore the 
surviving quantity in the last line of Equation (4.3b) is proportional to I/7 2 , and so the 
high-bombarding-energy dependence of a z (oo) is given by 



K — 

7C 



\ h 2 Ln( 7 )/7^ 



as was the case for the Goldhaber- Teller model. The different behaviors of a y (oo) and 
a z (oo) are reminiscent of the distinction between transverse and longitudinal "pulses" in 
the Fermi- Weizsacker- Williams method of virtual quanta ( 12 ]. Chapter 15). 

We see that in both the Goldhaber- Teller and Steinwedel- Jensen models, at sufficiently 
high bombarding energy, a z will become negligible compared to a y . The geometry of the 
collision also requires that a x = 0. In this situation, the argument of the Legendre poly- 
nomial in Equation (2.20) approaches |«y| 2 /|(ay) 2 | = 1, and the excitation probability of a 
state of specified n, i reduces to the simpler form 

' n ' e (2n)\\(2n + 2£ + l)\\ [lyl) ' 1 } 

Figure 3 illustrates the bombarding-energy dependences of a y and a z in the two models. 
The reaction involves 208 Pb projectiles and a 40 Ca target. The approach of a y to a constant 
limiting value, while a z strongly decreases, is evident. It is also clear from Figure 3, that all 
the cross-sections predicted by the Goldhaber- Teller and Steinwedel- Jensen models will be 
very similar, since all the cross-sections are determined by the two parameters a y and a z . 

We get a further simplification of Equation (4.4) if we restrict our attention to levels with 
a given total number of quanta N = 2n + £. Then we can deduce that 



VN fi/ ^ (2£ + l)(N + £ + 3) 



(T N-i n 

7^00 — 2 — 1* 



(4.5) 



V n-i-2 /12 (2£ + 5)(N~£) a N -i-2 /l2 

The second relation holds because the ratio is independent of a y , and therefore independent 
of b. According to Equation (2.22), if the ratio of excitation probabilities is independent of 
b, that ratio will also be the cross-section ratio. 
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In the particular case of 2 phonon levels, Equation (4.5) yields o%o/ &o,2 — ► 1/2- This is 
in agreement with the calculation of Bertulani and Baur |8| . 

Figure 4 shows the bombarding energy dependence of the excitation cross-section for 
levels with four or fewer GDR phonons in a 40 Ca target, when the projectile is 208 Pb. The 
calculation was done using the Goldhaber- Teller description of the GDR. As was shown in 
Section III, the predictions based on the Steinwedel- Jensen description would be similar. 

At bombarding energies per nucleon near 10 GeV, Figure 4 exhibits cross-section ratios 
consistent with Equation (4.5). At lower bombarding energies per nucleon, say below 1 GeV, 
the cross-section ratios are shown by Figure 4 to be quite different . Changing the bombarding 
energy has a significant effect on the ratio of transverse and longitudinal impulses received 
by the target. It follows from Equation (2.20) that this can affect the ratios of excitation 
probabilities of states of different angular momenta. Indeed, we see that as the bombarding 
energy per nucleon increases from 1 to 2 GeV, the N = 3 and N = 4 levels that are most 
strongly excited change from I = N to £ = N — 2. This behavior suggests that some 
interesting changes in the angular distribution of decay products might be observed as the 
bombarding energy moves through this region. 



APPENDIX A: PROOF OF THE POLYNOMIAL IDENTITY (2.14) 

We use the expansion 

(a • b) p = p\ V Jf f + \ PAk ■ b) 

e=p,p-2, P -A,... v ^ ' ^ ' 

in the power series for the exponential 
to get 



pi ^ — i p\ 



00 2/ -I- 1 

e-^EW £ (p _ Q „ (p + < + 1) „ ^-6) (AD 

p=0 l=p,p-2,p-A,... v ^ J \f ' ) 



£)(2£ + l)(a6)^(&.fi) J2 



(a 2 b 2 ) n 



n (2n)\\(2n + 2£ + 1)!!' 
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The Legendre polynomial can be decomposed by means of the spherical harmonic addition 
theorem 

i i 
{ab) e P t (a-b) = ^fp[(ab) e £ ^(«)^) = ^7 E (-l) m ^™(a)^(b), 

m=—£ m=—£ 

leading to 

" - <- £ ^ggM^ ± i-WM m 

n ,e=o K ' K ' m=-e 

which is Equation (2.14). 



APPENDIX B: THE NORMALIZATION FACTOR FOR 3-DIMENSIONAL 
HARMONIC OSCILLATOR STATES (EQUATION (2.15B)) 



We start with 



y*(4, c + 4)\o> 



^(4 + <)<|o> 



Define 



for which 



cT = 



c x + ic y 



V2 



c+ = 



[c+,4] = 1 , 

so that < 0|(c+) £ (c+) £ |0 > = t\. Thus 

< o| M(4,cf,c+)] + ^(4,^,4)10 > 



v/2 ' 
c + |0 > = 0, 



-1) £ /(2£+l)! 



An 



(Bl) 

(B2) 
(B3) 



* = (B4) 
4tt • 2^! 1 j 



.c / 2 + (4) 2 + (4) 2 - when ^ acts 



Now define the spherically symmetric operator H + = (c 
on a state it raises the number of quanta by two, without changing the state's rotational 
transformation properties. Using (2.12), it can be shown that 



[H,H- 



47V + 6, 



(B5) 



where N = c+c x + ctc y + c[c z is the number operator. Then 



H n (H + ry e (ct,4,ci)\0> = H n - 1 HH + (H + ) n - 1 y t e (4,<i,ct)\Q> 
H^H+HiH+y-^ict, c+c+)|0 > +(4(2n-2 + l) + Q)H n ~ 1 (H+y- 1 y e e (ct,4,ct)\0 > 



24 



If H is repeatedly commuted past factors of H + , the result is 

H»(H + ry t l (c+,<*,c+)\0> 
= [{4(2n - 2 + £) + 6} + {4(2n - 4 + £) + 6} + • • • {A£ + 6}] H n -\H + ) n - l y f e (ct , c+ c+)|0 > 
+ H n -\H+rHy*(4, c+ c+)|0 > (B6) 

But the last term in (B6) vanishes because Hyf(c£, c+, c+)|0 > would have to be a state 
with £ — 2 quanta, and angular momentum quantum nunbers (£,£), which is impossible. If 
we sum the arithmetic series in (B6), we get 

H n (H + ) n y" e (ct, c+,c 2 + )|0> = 2n(2n + 2£+l)H n -\H+r~ i y £ (ct,c+,4)\0> (B7) 

Iteration of (B7) gives 

^(^ + r^(4,<,4)|o> = (2n)!! ( ( 2 2 ^ 1 2 ) !, + m yi& 4, 4)|o > 

If this is combined with (B4), the result is 

<o\[(H+ry e (ct, c +, ct)} + (H+ry i (ct,4,ct)\o> = (2n)!!(2 " + 2f+1)!! , 

so that a normalized 3-dimensional harmonic oscillator eigenstate with angular momentum 
£, angular momentum ^-component £, and 2n + £ quanta can be written 

( ~ 1)W V / (2n)!!(2n 4 +2 lTTin ^ + (C * } ' + * > (B8) 

A factor (— l) n has been inserted to yield phases for the eigenstates consistent with those 
used by Brody and Moshinsky (fl?)]). 

Since the scalar products < \E^|\I/^ > (m = £, £ — 1 • • • — £) are independent of m, it 
follows that the normalization factor given in (B8) is also applicable if y\ is replaced by y^. 
This completes the justification of Eq.(2.15b). 

APPENDIX C: INTEGRATION ALONG THE HYPERBOLIC ORBITS 

In order to calculate the Fourier transforms of the expressions for F(t) and G(t) given in 
eqs. (3.18), we need the following integrals 



i(r'(t)} 3 loo [(r'(t)] 3 r'(t) 7-oo *"(*) 1 ' 
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It is convenient to express these integrals in terms of a parametric representation of the 



projectile orbit 



10|: 



y '(w) = a [cosh(» + e] (C2a) 



z'(w) = aVe 2 - 1 sinh(w) (C2b) 



a 



t( w ) = -[esinh(ty) + w] (C2c) 



v 



where w is a parameter (—00 < w < 00), v is the asymptotic projectile speed, and a and e 
are defined in III. 2. From (C2), it follows that 



r'{w) 


= a(l + e 


cosh(u>)) 


dt 


= —dw 

V 




y dt 


= dy' = 


a smh(w)dw 


■ 1 

z 


= dz = 


a\J e 2 — 1 cosh(w)d 



We can then see that 



00 



e^J^Ldt = r^ smHw)+wK sinh( W ) 



[(r'(t)] 3 av J _ 00 (1 + a cosh(w)) 2 



av 



e it{esinH W ) +W ]d f_h x + ecosh ( w ))-A dw 

7-oo dw V e / 



eat> y_ 00 aw v 7 



i£ — 1 / e i€[esinh{u,)+w] rfw (C3) 



eat> 



—00 



This last integral can be cast in a more convenient form by displacing the integration path 
from the real w axis to the line w + iw/2 (—00 < w < 00). We can do this without changing 
the value of the integral because the integrand has no poles between these two lines. Then 
(C3) becomes 

j^jS^dt = i£L± i / e W '-=$ '-)+v,+i*/2] dw 



00 



[(r'(t)] 3 * eav 



00 

00 



i ^ t2 V g " /2 I e -^ cosHw)+ ^ w dw 



eav 



-00 

2^- e~^ /2 / e- fecosh(M,) cos(^)dM; (C4) 
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To proceed, we define /(e) by 



m = 



it[esmh(w)+w] Sinh(w) ^ 

(1 + acosh(w)) 2 

t poo 

= 2i|e-^ /2 y e -«" osh (») cos(£w)dw 

Then a partial integration and some re- arrangement show that 

cosh(«,)+» e + cosh(w) ^ 
(1 + ecosh(w)) 2 

POO 

-2e^ n/2 / e^ ecosHw) cos(w)dw 
Jo 



d_ 

de 



d_ 

de 

2£e" 



-Z-k/2 



3 — £e cosh(wj) 



cos(w) cosh(w) dw 



Thus 



j — t 



[(r'(i)P 



1 [°° e it[e S in Hw)+W ] e + COsh( W ) ^ 

(1 + ecosh(w)) 2 



av 



= ^ e -^/2 /" e-^" 0811 ^ cos(«;) cosh(w)rfw (C5) 
av J 



iwt V'jt) fa = f°° ^[esmh(w)+w] sinh(w) 



This same shift of the integration path leads to 
and 



1 + e cosh(w 



-dw 



^ e-*""*™ cosh( W ) C ° s( ^ ) + e ^ Hw) ^ w) dw(C6) 
Jo 



1 + (e sinh(w)) 2 



e i£[e S inh(w)+w] COsll(w) ^ 



1 + e cosh(w) 



= 2V^le-^ r e-^^H sinh( W ) = +/ ^\^H^ 

Jo 1 + (esmh(w)) 2 

If the four integrals (C4, C5, C6, C7) are used in (3.18) and (3.8b), the result is the set 
of expressions (3.19) for a y and a z . 
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Figure 1 

FIG. 1: The vector R connects the neutron mass center to the proton mass center, locates the 
i th proton relative to the nuclear mass center, whereas p'j locates it relative to the proton mass 
center. / is the ratio iVr/M-T- 




Figure 2 



FIG. 2: Comparison of results of relativistic and non-relativistic calculations of the cross-section for 
excitation of the one-phonon GDR level in 40 Ca, due to Coulomb excitation by 208 Pb projectiles. 
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FIG. 3: Comparison of calculations of (a y ,a z ) using the Goldhaber-Teller and Steinwedel-Jensen 
models of the GDR. a) b = 12 fm; b) b = 45 fm. 
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FIG. 4: Excitation cross-sections for various (n,£) levels of the GDR in 40 Ca, due to Coulomb 
excitation by 208 Pb projectiles. Levels with the same number of phonons are indicated by the 
same type of line. These calculations are done using the Goldhaber- Teller description of the GDR. 



